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C\i ■ Abstract 

Q ■ The QCD description of exclusive decays of the A& baryon involves hadronic matrix elements 

of non-local light ray operators, the light-cone distribution amplitudes. We introduce the 
complete set of three-quark distribution amplitudes and calculate the renormalization scale 
dependence for the leading twist. At leading order in the strong coupling the evolution is 
driven by pairwise two-quark interactions: heavy-light involving Sudakov logarithms as in 
the B-meson case, and light-light as in light mesons. We solve the evolution equation and 
show that its main effect is to generate a radiative tail extending to high energies. Finally, we 
present simple models for the distribution amplitudes based on QCD sum rules, and study 
the effect of the evolution on these models. 
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1 . Heavy baryons containing a 6-quark will be copiously produced at the LHC and their 
weak decays may provide important clues on flavour-changing currents beyond the Standard 
Model. A particular advantage of A& baryons over B mesons is their spin. Their polarization 
facilitates the study of spin correlations, providing valuable information on the chirality of 
the short-distance transition. This can be exploited for example at LHCb by studying rare 
radiative decays such as A b — > A7 [1]. 

The theory of 6-baryon decays into light hadrons is, however, more complicated and was 
receiving less attention compared to the B-meson decays. In particular, we are not aware of 
any dedicated study of the heavy-baryon distribution amplitudes (DAs) that are the primary 
non-perturbative objects required for calculating heavy-baryon decays into light particles 
based on the heavy quark expansion, see e.g. [2], or using sum rules of the type proposed 
in [3-5]. The only existing models of heavy baryon DAs [6,7] are motivated by quark models 
and are not consistent with QCD constraints. In this letter we give, for the first time, the 
complete classification of three-quark DAs of the A& baryon in QCD in the heavy quark limit 
and discuss some of their main features. In particular we derive a renormalization-group 
equation that governs the scale-dependence of the leading-twist DA and study its solution. 
Simple models of the DAs are suggested, and their parameters are fixed based on estimates 
of the first few moments using QCD sum rules. 

2. The A{> distribution amplitudes can be defined as matrix elements of non-local light- 
ray operators built of an effective heavy quark and two light quarks following, on one hand, 
the similar construction for B-mesons [9-14] and, on the other hand, the QCD description 
of nucleon DAs [15, 16]: 

e abc (0\ (u a (hn)C l5 iid b (t 2 n)) h c v (0)\A(v)) = /f * 2 (M 2 )A(i;) , 
e a6c (0| (u a (t in )C l5 d b (t 2 n)) K(0)\Hv)) = f^Kh, t 2 )A(v) , 
(0| (u a (t in )C l5 ia nn d b (t 2 n)) feS(0)|A(u)> = 2fg ) *S{t 1 ,t 2 )A(v) , 



abc 



abc 



(0| ( M a (tm)C75^(t 2 n)) h c v (0)\A(v)) = /f* 4 (<i,t2)A(T;) . (1) 



Here the subscript 2, 3, 4 refers to the twist of the diquark operator, C is the charge con- 
jugation matrix, A(v) is the Dirac spinor, rfA(v) = A(v). The non-relativistic normalisation 
is assumed, AA = 1. Further, n M and n M are light-like vectors which we choose such that 
^ = Ov + ^m)/ 2 > v -n = l,n-n = 2. 

The DAs \J/ 2 , ^35 ^4 are symmetric under the interchange of the coordinates of the two 
light quarks; ^3 is antisymmetric. The couplings /| are given by the matrix elements of 
the local operators 

e a6c (0| (u a (0)C l5 d b (0)) K(0)\Hv)) = f^Aiv) , 

^abc 



e 



(0| (u a (0)C 75 K(0)) fcS(0)|A(v)) = ffA(v) . (2) 

and are used in ([T]) for convenience. 



Alternatively, one can define 

e abc (0\ (« B (tin)C75^(t 2 n)) ^(0)|A(i;)) = f^fihMH'") 



abc 



(0| (u a (t in )C l5 $ild b (hn)) ^(0)|A(v)> = /i 1) *J + (*i,t 2 )A(7;) 



(3) 



In contrast to the above, the DAs in (j3j) do not have any definite symmetry; however, the 
isospin zero condition for the diquark implies that 



It follows that 



% + (h,t 2 )=Vt~(t 2 ,t 1 ). 



*3(*1,*2) = T ^"(*1,*2)+*3"(*2,*1 



(4) 
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*r(*i,*2)-*r(fe,*i) 



(5) 



correspond to the symmetric (antisymmetric) parts of ^3 , respectively. 
Going over to the momentum space, we define 



so that 



00 /-00 

dwi / du 2 e" 
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-itltdl— it2UJ2 



1J)(U\,(JJ2, 



00 />1 

ujduj j due- iuj{tlU+it2ii) ip(uj,u),(Q) 
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ip(u,u) — ip(uuj,vuj) (7) 

where u = 1 — u. In the first representation u>i and u;2 are the energies of the u- and <i-quark, 
respectively, and in the second one uj = uj\ + u 2 is the total energy carried by light quarks (in 
the heavy-quark rest frame) whereas the dimensionless variable u corresponds to the energy 
fraction carried by the w-quark, i.e. uo\ = uuo and u 2 = uuj. These two representations are 
fully equivalent and can be convenient in different contexts. 



(2) 



„(1) 



A non-relativistic constituent quark picture of the A b suggests that /^ — f& at low 
scales of order 1 GeV, and this expectation is supported by numerous QCD sum rule cal- 
culations [19-22]. In fact, the difference between the two couplings is only obtained at the 
level of NLO perturbative corrections to the sum rules [21,22] and it is numerically small. 

The anomalous dimensions of the operators in (J2J) 



din /1 
are known to NLO 



"7 



(0 



$^a(/i) 



4 4) ; 



a(Ai) = af 0i)/(47r) 



(8) 



7i 



(i) 



7i 



(2) 



7$" 



(2) 
72 



796 - 16C(2) - 40n 7 
322 - 16C(2) - 20n ; 



(9) 
(10) 



§72 quoted here are in the naive 75 scheme. 



Thus, the scale dependence of the couplings is given by 



r(0 



«/ 



/A V) = WM 






TfV/3o 



a s (n ) -a s (/i)7j ; / 7^ 



W 



W 



47T 



/3 



v7i 



W 




(11) 



where do(//) /d In /i = -/3 Q a(ij,) 2 -l3 1 a(fi) 3 + - • • with /3b = 2(11 -2/3n/), A = 4(51-19/371/). 

For the numerical value of the couplings we quote the result of the NLO QCD sum rule 

analysis in Ref. [22]: 



f (2) 



./, 



(i) 



0.030 ± 0.005 GeV 3 



(12) 



at the renormalization scale /1 = 1 GeV. Note that these couplings cannot coincide at all 
scales since the corresponding operators have different anomalous dimensions. 

Similarly to the B-meson case [10,14,17] QCD equations of motion can be used to derive 
exact relations between the three-quark DAs in ([1]) and the four-particle DAs involving an 
extra gluon field strength tensor. The corresponding analysis will be presented elsewhere. 

3. The DAs in ([TJ are scale dependent. The leading-order (LO) evolution equation for 
the leading- twist DA fai^i, uj 2 ; /i) can be derived following the usual procedure by identifying 
the ultraviolet singularities of one-gluon-exchange diagrams. 

The result can be expressed in terms of the two-particle kernels familiar from the evolution 
equations of the .B-meson and 7r-meson DAs. We obtain 



tlj-lfa(Wl,W2',V) 




du[-f {u/ 1} Ul]fi)lff2{u' 1} U 2 ]fi) 



+ 



d^2l 04,^2; V)i>2{Vl,W2,fJ>) 



- / dvV(u,v)^ 2 (vuJ,vuj;n) + -i)2(^i,^2^) 



(13) 



where in the last line u = u>i + u>2 and u = u>x/(u>i + 002)] the last term in the curly 
brackets, \ip2, is a result of the subtraction of the one-loop renormalization of the coupling 
f K according to Eqs. (JHJ) and ffTDl . 



The first two convolution integrals in Eq. ( 1T31) are associated with heavy-light dynamics: 
each of them involves just one of the light quarks. Indeed, the kernel 7 LN (u/, u>; //) coincides 
with the one controlling the evolution of the B-meson distribution amplitude, the Lange- 
Neubert anomalous dimension [11] 



7 L V,^;/i) = 
r LN (u/,u;) = 



In ) 6(u - u') - r LN (o; / , u) 

u> 4 



uo 6(u' - w) 6(uu - u') 



CO' LU' 



uJ 



LU — LO' 



where 



duo' 



7(««/,w) /(wO 



du'~i(u\uj) f(uj') -/(w) 



(14) 



(15) 



In turn, the last convolution integral in Eq. (TT3l) describes the interaction between the 
light quarks. V(u,v) is the celebrated ER-BL kernel [8]: 



V{u,v) 



u - v) + - H 0(v - u) 



1 — V \ u — V ) V \ v — u 

where the "+" subtraction is defined as 



(16) 



[V(u,v)] + = V(u,v)-5(u-v) I dtV(t,v). (17) 

Jo 

Note that in Eq. (fl3|) we retain the dependence on the number of colors N c in the prefactor 
although the whole construction only makes sense for N c = 3. 

4. For small evolution ranges, ln(/x//xo) ^ 1, it is sufficient to interpret the derivative on 
the l.h.s. of f|T3|) as a finite difference [ip 2 {^i, ^2\ A 4 ) — ^(^i, ^^Lpo)]/ ^(aVMo) an d substitute 
the initial condition ^2(^1, ^2; /-to) f° r ^2(^1, ^2; A 4 ) on the r.h.scl Obviously, this corresponds 
to taking into account one-loop renormalisation only, neglecting the resummation of poten- 
tially large logarithms. As we shall see below (Figures [2] and |3]) this single-evolution-step 
(one-loop) approximation is quite good in practice, e.g. for fj, = 1 GeV and \x ~ m&/2. 

In order to go beyond the one-loop approximation, one possibility is to integrate the 
evolution equation (fl3l) numerically. We have taken another, semi-analytic, approach which 
has an advantage that it allows one to understand the structure of the solution. To this end 
we first remove the ln(/i) term on the r.h.s, which is related to the cusp anomalous dimension, 
by defining: 

^2(^1, w 2 ;^) = #*>i, w 2 ; A*) f — 2^) , (18) 

where 

f 1 dm CFOLsim) 
g(a s (fjL))= T cusp (a s {m)), T CUEp (a s {m)) = 1 . (19) 

7/10 m An 

Substituting (|T8"|) in (IT51) yields an evolution equation for <f>(u)i, uj 2 ; A 4 )- The next step is to go 
over to the moments space: 

POO 

.JV-1 / j, . , .Af-1 



<p(N,M-,(i)= dujiLu?- 1 du 2 u 2 M - l 4){uj 1 ,uj 2 \n). (20) 

Jo Jo 

This leads to factorization of the Lange-Neubert terms since 

duu N - 1 T LS {u , 1 u) = (^ , ) iV - 1 f LN (iV); f LN (iV) = -tt(JV) - *(-iV) - 2 7B . (21) 



'0 
The ER-BL term does not factoriz 



f*(JV,M;/x) = 



<it> / <iw / du;u;(uu;) (ucj) ' ( - r J V(«, v) (J)(vuj,vuj; jj) 



a jo jo ^uuJ 



<KN,M;n) 

(22) 



^Note that the scale [i appearing explicitly in (TTJ} and the scale of the strong coupling must be the same. 
"We use the shorthand notation u = 1 — u, v = 1 — v. 



calling for some approximation. A simple one is obtained by substituting the initial condition 
4>(vu), vu; /i ) for (j)(vu,vu; fi) in both the numerator and denominator of (1221) . While this 
can be a starting point for an iterative procedure, we find that in practice such iteration is 
not necessary owing to the smallness of the ER-BL term. 

With this assumption, the r.h.s. of the evolution equation for <p(N, M; /i) factorizes leading 
to exponentiation of the kernels: 

f f dm 1 

(f>{N, M; fi) = <f>(N, M; /i ) exp / E{N, M, m) \ (23) 



where the exponent is given by: 

E(N,M,m) = g(a s (m)) + ' 2nJm) 



l + r LN [N-g(a s (m))/2 



3n 



+ f LN [M-g(a s (m))/2 ) +t^(N,M;m) 



(24) 



resumming double- as well as single-log terms to all orders. 

Further simplification in evaluating Ty(N, M; fj) is achieved by replacing the ER-BL 
kernel V(u, v) in (|22l) by the expansion as a sum of products of Gegenbauer polynomials (see 
e.g. [24]): ' 

V(u,v) = -u(l -u)JT 2 i 2 ^ 3 l^ nC 3 J\2u - l)CT{2v - 1) , (25) 

^ ln+ l)(n + 2) 

where 7 n = 1 — 2/[(n + l)(n + 2)] + ^YTm=2 ^/ rn is the leading-order anomalous dimension. 
In the numerical evaluation presented below we truncate the sum in f )25p at the leading 
non-trivial term, n = 2. The impact of this truncation proves to be small, at least for the 
models of the DA that we consider. 

Finally, the answer in the momentum (energy) space is restored by a double inverse- 
Mellin transform which is done numerically. The main effect of the evolution is to generate a 
'radiative tail' of the DA that falls off as 1d.(u>i/ju)/u>i or Id^u^/j^)/^ & t large energies, which 
is the same effect that the evolution has on the B-meson DA, see [11-14]. 

5. Realistic models for the DAs can be obtained using QCD sum rules for the correlation 
functions involving the non-local light-ray operators in ([T]) and a suitable local current. We 
define 

J(x) = e abc (d a (x)P + -f 5 C T u\x)) h c v (x) (26) 

where P + = (1 + i>)/2 and consider, for the leading twist, the correlation function 

d 4 xe- tEvx e abc (0\(u a (t 1 n)C-f 5 '4d b (t 2 n))h c v (0)J(x)\0) = 

= P+ uduj due- luj{utl+ut2) U 2 (iu,u;E) (27) 

Jo Jo 



and similarly for the other structures. The general form of the sum rule is then 

i|/f \ 2 ^ R (co,u)e~^ = M[U 2 }(u;,u;r,s ), (28) 

where M[U2](uj,u;t,Sq) is the Borel-transformed continuum-subtracted invariant function 
n 2 (co', u; E); r is the Borel parameter which we take to be in the interval 0.4 < r < 0.8 GeV 
and So = 1.2 GeV is the continuum threshold (interval of duality); A = m^ b —nit, ~ 0.8 GeV. 
Taking into account only the leading-order perturbative contribution to the sum rule, 
one obtains 



^ATVim [ ° dse- s/T (s - u/2) 

2 Jw/2 



with 



i/j 2 {u,u 

'»/ 

^ A (u, u) = 5AT 1 f ° ds e- s/T (s - uo/2f , 

Jw/2 

V»5(w, u) = —M^uo / ds e- s/T (s - u/2) 2 , 

4 Jw/2 

^(w, u) = ^Af~ 1 iu(2u - 1) / ° ds e~ s/T (s - u/2) 2 , (29) 

4 Ja;/2 



Af = I °dss 5 e- s/T . (30) 



To this accuracy the coupling is equal to |/a| 2 = e A//r A/"/(207r 4 ). All DAs have in this 
approximation the support property < uj < 2so and are normalized such that 

" 2s o r 1 „ / >2s o /•! _ /-2so z-i _ 

cjcL; / dutl>2(u,u) = / wdo; / du^(u>,u) — / u;gL> / duip^u^u) — 1 . 
Jo ./o ./o ./o ./o 

(31) 

Note that the leading-twist DA ip 2 (u,u) = ip 2 (u)i, ^2) vanishes when either one of the light- 
quark energies goes to zero: lj 2 u(1 — u) — uj\uj 2 . This property is model-independent and 
consistent with the evolution equation in ffT3|) . 

The limit r — > 00 is known as the approximation of local duality. In this case one obtains, 

for example 

45 

-02 D (^i, w 2 ) = ^-6^i^2(2s - Wi - lu 2 ) 2 9(2s - wi - cj 2 ) (32) 

5S 

and similarly for other twists. With decreasing Borel parameter the DA becomes tilted 
towards smaller momenta since contributions of large uj —>■ 2so are more strongly affected by 
the additional (exponential) suppression factor. 

As explained in [9,12] in order to evaluate the non-perturbative contributions to the sum 
rule one is forced to use the non-local condensates. We use the general parametrisation [25,26] 

poo 

(q(x)q(O)} = (qq) / due^^f(u) (33) 

Jo 
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where (qq) ~ —(240 MeV) 3 is the quark condensate and f{y) is the model function [12,27] 

/(*0 = ^ \y - a e- x/v ; a-3 = 4; (34) 



no -2i 



m 



2 ' 

o 



m,Q ~ 0.8 GeV 2 is the standard notation for the ratio of the mixed quark-gluon and quark 
condensates, and A ~ (400 MeV) 2 is the correlation length. 
Using this model, we obtain the sum rule: 



il/fWV.tOe-^ 



167T 4 



UJ 2 UU 



dse- s/T (s-uj/2) 



u/2 



(uu) u a _ 2 sin got 



iir 2 u 7r(3 — a) 



(s -iu/2-K) 3 - a e- So/T 



+ - / ds e~ s/T (s - uj/2 - k) 



3—a 



(dd) u _ a _ 2 sina7r 
87r 2 u 7r(3 — a) 



(s -u/2-R) 3 - a e- So/T 



1 /" s ° 

+ - / dse- s/T (s-uj/2-R)) 



— \ \3— a 



- (mm) {dd)T 2 f(2Tuuj)f{2Tuu)e- ul/{ - 2T) 



(35) 



where Heaviside functions of the difference between the upper and the lower limits of inte- 
gration are implied, and where we used a shorthand notation 



A 



K 



2uu' 



A 

2uu 



(36) 



From the vast experience of QCD sum rule calculations of the pion DA (see e.g. [25,28-30]) 
it is known, however, that the QCD sum rules cannot give the functional form of the DAs 
but rather have to be used to constrain certain momentum fraction integrals (the moments). 
Furthermore, obtaining a meaningful error estimate is especially difficult in the present case 
because there is not enough experience in using the concept of non-local condensates in 
baryon sum rules. In order to be on the conservative side we adopt the following procedure. 
As well known, QCD sum rules can be written for different interpolating currents. Our choice 
in (|2"B"|) corresponds to the constituent type sum rule, in the terminology of Ref. [22], and 
has the advantage that the corresponding sum rule (135]) has several terms. Replacing the 
projector P + in (|26|) by the unity matrix or by js one obtains two other currents, called J\ and 
J 2 in [22]. Obviously J = {J\ + J-i)j2. The corresponding sum rules pick up contributions 
of even and odd dimension in (I3"5"|) respectively, i.e. perturbation theory and the quartic 
condensate for J 2 and the quark condensate for J\. We take the sum rule in (133|) for our 
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Figure 1: QCD model for the leading twist DA of the A(, baryon defined in Eq. ([38 



central values and use the spread of the results using J\ and J 2 (using the central values of 
the parameters) as an error estimate: 



2s o /-i _ 

ujduo I duip2(u,u) = 1, 
n Jo 

dm I duil)2{oJ-,u) = 1.7 ±0.7. 
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-,3/2, 



2s /-l 
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2s rl 



0.5 
) 



n 



^3/2 



+0.7 



du / duC 6 2 '\2u- 1) ^ 2 (cj,m) = 0.6±i 



(37) 



The error bands given in f|3T|) should be regarded as most conservative: using the sum rule 
(131)1) alone and varying the parameters in a reasonable range yields much smaller variations. 
In particular the value of the first integral in the second line in (1371) is very stable with 
respect to variations of the Borel parameter. 

Note that the ratio of the integrals with and without the u factor is different for the 
Gegenbauer moment in the energy fraction as compared to the unit weight. This implies 
that the w-dependence of these components is different. Taking into account the expected 
low-energy behaviour ~ iO\Ui and the sum rule moments of Eq. (1371) . we propose a simple 
model (see Fig. 1) for the leading-twist DA at the low scale of /i = 1 GeV: 



^2(^5 u ) — oj 2 u{\ — u) 



K-^ 0+a2 c^ 2 (2u-l)\e 







(38) 



>+0.250 



with e = 200+H MeV ; £ i = 650ijjj;g MeV and a 2 = 0.333i _.,.,.,. 

In the calculations of A^ decays into light quarks using QCD factorisation one expects 
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Figure 2: QCD model for the leading-twist DA of the A& baryon defined in Eq. (|38|) at the scale 
of 1 GeV (solid curve) and after the evolution to [i = 2.5 GeV (dash-dotted curve) as a function 
of uj = uj\ + u>2 for two values of the light quark momentum fraction u = 0.5 and u = 0.125. The 
result of a single-step evolution to \i = 2.5 GeV, which includes the ~ 0(a s ) correction only, is 
shown by dashes for comparison. 



that integrals involving negative powers of the quark momenta will contribute, for example: 

6>(A UV - u\ - uj 2 ) 



/•Auv /"l ^ u _ /"co rco 

A g (//, Auv) = / du; / — ^2(1^, u; /j,) = / dwi/ cL> 2 
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A d (/x,A uv )= / cL> / duip 2 (u,u;iJ,) 



00 />oo 



. 6>(A UV - wi - u 2 ) , 
dut / du 2 ■ ip 2 (ui , v 2 ; (A) , 

JO U)\-\- L0 2 

(39) 

where an additional energy cutoff u < Auv is introduced in the definition of the moments. 
This guarantees that the moments are finite in presence of a radiative tail - - such a tail 
will be generated by evolution to a higher scale, even if it is not introduced at the low scale 
in a given initial-condition model. We recall that in the B-meson case the radiative tail 
~ ln(u) / fj) / u renders such an energy cutoff necessary for any positive moment, while the 
first negative moments analogous to ( 1391) would be finite in its absence. The two-dimensional 
integration in the A& case implies greater sensitivity of the cutoff, as even the first inverse 
moment A q diverges in its absence. 

For the above model without an energy cutoff one obtains: 



A,(l GeV) = - + ^ ~ 5.5+gJ GeV" 1 ; 

En £1 



A d (lGeV) = — ~ 1.7 ±0.7 GeV- 1 . (40) 
3e 



Since £\ 3> Eq, the contribution of the Gegenbauer correction to the first integral is small 
so that a 100% uncertainty in a 2 does not play a significant role. The effect of an energy 
cutoff Auv °f order 2-3 GeV is already small. The central values of A d and A q for Auv = 2.5 
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Figure 3: The u dependence of the DA for fixed uo = 0.5 GeV (near the peak in Figure [2]) and 
lo = 1.0 GeV (crossing over to the tail region). The curves are as explained in Fig. 2. Note that 
the effect of evolution to higher \x is to decrease the DA for any u in the former case and increase 
it in the latter. 

GeV are summarized in Tabled], where we also present the renormalization scale dependence 
of these moments. The effect of evolution from li = llq = 1 GeV to Li — 2.5 GeV on the 
functional form of the DA is illustrated in Figures 2 and 3. 






/f [GeV 3 ] 


Ad [GeV- 1 ] 


A, [GeV- 1 ] 


li = 1 GeV 


0.0300 


1.66 


5.38 


Li = 1.5 GeV 


0.0314 


1.52 


4.94 


Li = 2.0 GeV 


0.0322 


1.40 


4.61 


Li = 2.5 GeV 


0.0329 


1.31 


4.34 



Table 1: The decay constant /^ of Eq. ([2]) and the typical integrals A^, A q of Eq. (|39|) at different 
renormalization scales li. The moments are all computed at a fixed energy cutoff Auv = 2.5 GeV. 
These numbers correspond to the central values of the model described above — the theoretical 
uncertainty is as quoted in Eqs. (fT2j) and (|4T)|) . 

A similar analysis can be done for the twist-three DAs. Without going into details we 
present the simplest models that are consistent with the QCD sum rule constraints: 



oo 



ip s 3 {uj,u) =—^e 



-ui/e 3 



2e| 



^a(u;,u) 



CO 



2eV 



(41) 



with e 3 = 230 MeV. 
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